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Abstract. We study the evolution of the one dimensional periodic cubic Schrodinger 
equation (NLS) with bounded variation data. For the linear evolution, it is known that 
for irrational times the solution is a continuous, nowhere differentiable fractal-like curve. 
For rational times the solution is a linear combination of finitely many translates of the 
initial data. Such a dichotomy was fist observed by Talbot in an optical experiment 
performed in 1836, 20 . In this paper we prove that a similar phenomenon occurs in the 
case of the NLS equation. 



1. Introduction 



In a 1836 optical experiment Talbot [20] observed white light passing through a diffraction 
grating. He looked at the images that were produced with the help of magnifying lens and 
noticed a sharp focused pattern with a certain periodicity depending on the distance. 
Rayleigh [18J calculated the Talbot distance as ci = ^ where a is the spacing of the grating 
and A is the wavelength of the incoming light. 

Berry with his collaborators (see, e.g., [H 121 [Sj H]) studied the Talbot effect in a series 
of papers. In particular, in Berry and Klein used the linear Schrodinger evolution to 
model the Talbot effect. They showed that at rational multiples of the Talbot distance 
overlapping copies of the grating pattern reappear, while at irrational times the images 
have a fractal nowhere differentiable profile. Also in [Ij, Berry conjectured that for the 
n— dimensional linear Schrodinger equation confined in a box the imaginary part Qu{x,t), 
the real part ^u{x, t) and the density \u{x, t) P of the solution is a fractal set with dimension 
D = n + for most irrational times. He also observed that in the one dimensional case 
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there are space slices whose time fractal dimension is | and there are diagonal slices with 
dimension |. 

The first mathematically rigorous work in this area appears to be due to Oskolkov. In 
|17j . he studied a large class of linear dispersive equations with bounded variation initial 
data. In the case of the linear Schrodinger equation, he proved that at rational times the 
solution is a continuous function of x and at rational times it is a bounded function with 
at most countably many discontinuities. For the exact statement see Theorem [3] below. 

The idea that the profile of linear dispersive equations depend on the algebraic properties 
of time have been further exploited in the papers of Kapitanski-Rodniaski [H], Rodnianski 
|19j . and Taylor [21]. In [H], the authors show that the solution to the linear Schrodinger 
equation has better regularity properties (measured in Besov spaces) at irrational than 
rational times. It is important to note that this subtle effect can not be observed in the 
scale of Sobolev spaces since the linear propagator is a unitary operator in Sobolev spaces. 
In [19], using the result in [T3], Rodnianski partially justified Berry's conjecture in one 
dimension, see Section [21 

In [21], Taylor independently obtained Berry's quantization result and also extended it 
to higher dimensional spheres and tori [22]. In particular he proved that at rational times 
the solution is a linear combination of finitely many translates of the initial data with the 
coefficients being Gauss sums. He further showed that some classical identities for Gauss 
sums can be obtained by an analysis of the linear Schrodinger evolution. This shouldn't 
come as a surprise since number theoretic connections of Talbot effect had already been 
observed in the works of Oskolkov, Kapitanski and Rodnianski. In [21], Taylor also noted 
that the quantization implies the boundedness of the multiplier e**^ for rational values 
of It is known that, [21], the propagator is unbounded in for p 2 and ^ irrational. 
This can be considered as another manifestation of the Talbot effect. 

More recently Olver |16] , and Chen and Olver [HI E] provided numerical simulations of 
the Talbot effect for a large class of dispersive equations, both linear and nonlinear. In the 
case of polynomial dispersion, they numerically confirmed the rational/irrational dichotomy 
discussed above. This behavior persists for both integrable and nonintegrable systems. 
An important question that the authors raised is the appearance of such phenomena in 
the case of nonpolynomial dispersion relations. The numerics demonstrate that the large 
wave number asymptotics of the dispersion relation plays the dominant role governing the 
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qualitative features of the solutions. We should also note that in [2^ the Talbot effect was 
observed experimentally in a nonlinear setting. 

In this note we investigate the Talbot effect for cubic nonlinear Schrodinger equation 
(NLS) with periodic boundary conditions. Our goal is to extend Oskolkov's and Rod- 
nianski's results for bounded variation data to the NLS evolution, and provide rigorous 
confirmation of some numerical observations in [161 El IH]- In particular we prove that for 
a large class of rough data the solution of the NLS equation, for almost all times, is a 
continuous but fractal-like curve with upper Minkowski dimension|j |. Our main result is 
the following theorem: 

Theorem 1. Consider the nonlinear Schrodinger equation on the torus: 

iut + u^x + \u\^u = 0, t e M, a; e T = M/27rZ, 
0) = g{x). 

Assuming that g is of hounded variation, we have 

i) u{x,t) is a continuous function of x if ^ is an irrational number. For rational values 
of the solution is a bounded function with at most countably many discontinuities. 
Moreover, if g is also continuous then u G C^C^. 

ii) If in addition g ^ Ue>o-^^'''^' then for almost all times either the real part or the 
imaginary part of the graph ofu{-,t) has upper Minkowski dimension |. 

We note that the simulations in [161 El H] were performed in the case when g is a step 
function, and that Theorem [1] applies in that particular case. 

To prove Theorem[T]we first obtain a smoothing result for NLS stating that the nonlinear 
Duhamel part of the evolution is smoother than the linear part by almost half a derivative. 
For bounded variation data, this implies that the nonlinear part is in which immedi- 
ately yields the upper bound on the dimension of the curve. The lower bound is obtained 
by combining our smoothing estimate with Rodnianski's result in [19], and an observation 



^Upper Minkowski (also known as fractal) dimension, dim(i?), of a bounded set E is given by 



lim sup ■ , 
.->o' log(i) 



where M{E, e) is the minimum number of e-balls required to cover E. 
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from [Hj connecting smoothness and geometric dimension. We remark that the first part 
of Theorem [1] was observed in [12] in the case of KdV equation. 

2. Discussion of earlier results and the proof of Theorem [T] 
First of all recall that for s > 0, H^{T) is defined as a subspace of via the norm 

WfWH^m := I^Jky^lfm^ 

where {k) := (1 + k'^Y^'^ and /(A;) = j^^ f{x)e~'^^^dx are the Fourier coefficients of /. We 
note that the local and global wellposedness of NLS for data for s > was obtained by 
Bourgain 

We start by formally decomposing the solution as 

where P = \\g\\'^/T^- Here is the nonresonant part of the nonlinear Duhamel term of the 
solution. The proof of Theorem [1] relies on earlier results by Oskolkov and by Rodnianski 
on the linear part e*^^"^"'"'"^)*^' and the following smoothing theorem involving M . 

Theorem 2. Fix s > 0. Assume that g G H^{T). Then for any a < min(2s, 1/2), we have 

Af{x,t) e Cf^^H^t!^. 

We note that this type of smoothing does not hold for s < 0. For this and (local) smoothing 
results in J^i^ spaces, see [lO]. We will prove Theorem [2] in Section |3] below. Note that if 
g is of bounded variation then g G ne>o H^~^, and hence 

Ar(x,t)Gf|Q6M^.'a- 

e>0 

In particular, M^x.t) G C°C°. This and the following theorenj^ of Oskolkov [T7] conclude 
the proof of part i) of Theorem [H 

Theorem 3. [T7] Let g he of hounded variation, then e^^^'^'^g is a continuous function of x if 
^ is an irrational numher. For rational values of it is a hounded function with at most 
countahly many discontinuities. Moreover, if g is also continuous then e**^^"'(7 G C^C^. 

^In fact Theorem [3] is a special case of a Theorem of Oskolkov which asserts the statement above for 
one dimensional dispersive equations with polynomial dispersion relation. 
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Remark 1. In the case when the initial data has jump discontinuities, by the quantization 
result in [21 EH [16], the linear solution has jump discontinuities when ^ is a rational 
number, since it is a linear combination of finitely many translates of the initial data. 
Noting that is a continuous function of x, we have the same conclusion for the NLS 
evolution. 

The second part of the theorem will rely on a result of Rodnianski [19j which in turn 
relies on results obtained by Kapitanski and Rodnianski [14J. To state this result we need 
to define the Besov space -Bp^oo the norm: 

||/||B|^:=sup2^^||P,./|U„ 

where Pj is a Littlewood-Paley projection on to the frequencies ~ 2^ . We should note that 
C"(T) coincides with B^^, see, e.g., 

Theorem 4. [19] Assume that g is of bounded variation and that g ^ [J^^qH^^^, then for 
almost all irrational 

e^^^^^g^ljBlt:, 

e>0 e>0 

We note that for each a, H"" G Sf „„, and that, for < a < i, C C". Therefore 

by Theorem [H we have for each t, 

A/'(-,^)en^i- 
Af{;t)ef]c'^-\ 

This implies that Theorem H] is valid for the nonlinear solution u. 

It is a well-known result that if / : T — )■ R is in C°, then the graph of / has upper 
Minkowski dimension D < 2 — a. Therefore, the graphs of 3fJ(-u) and Q'(m) have dimension 
at most |. We note that this upper bound wouldn't have followed if we had less than half 
a derivative gain in Theorem [21 

The lower bound follows from the claim of Theorem [H for the nonlinear solution u and 
the following theorem of Deliu and Jawerth [11] . 
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Theorem 5. |TT] The graph of a continuous function / : T — )■ M has upper Minkowski 
dimension D > 2 — s provided that f ^ IJe>o ^i^- 

3. Proof of Theorem [2] 

The theorem follows from a simple X^''' space estimate that appears to be new. 
With the change of variable u{x,t) — )■ u{x,t)e^^'^, where P = \\g\\2/TT, "we obtain the 
equation 

iut + Uxx + \u\'^u — Pu = 0, t G M, X G T, 
with initial data in g E iJ*, s > 0. We want to prove that for any a < min(2s, we have 

„ _ p^dxxt rz TTS+a 

Note the following identity which follows from Plancherel's theorem: 



(k) = ^ u{ki)u{k2)u{k - ki + ^2) 



ki,k2 



— \\u\\lu{k) — \u{k)\'^u{k) + u{ki)u{k2)u{k — ki + k2) 

=: Pu{k) + f^u){k) + Riu){k), 

2-K 



where u{k) = J^^ u{x)e ^'^^dx are the Fourier coefficients. Using this in the Duhamel's 
formula, we have 

u{t) = e**^^^^ + i [ e^(*-")^^^(p(M) + R{u))dT. 







We note that 



:i) mu)\\h^^^ = Jj2\m\'{k)''^"' <\\urHs 



for < a < 2s. 
Using ([T]), we have 



\\u{t)-e''^^^g\\Hs+a< [ \\u{T)\\lsdr + [ e'^'-^^^- R{u)dT 

Jo Jo 

To bound the integral involving R{u), we work with the X'''^ space [5l |6]: 

llMllx^.b = \\u{r, k){ky{T - k"^) 
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with b > 1/2. We also define the restricted norm 

ll-ull^s.ft = inf IImII 

s u=u on [— <5,<5] 

We will use the embedding X^'^ C CfH^ for 6 > | and the following inequality from |5]. 

For any s G M, 5 < 1, and 6 > |, we have 

ft 



I 5-, \\F\\ ^s,b-i. 



Thus for < t < (5, we have 
ft 



J{t-T)dx 



R{u)dT 









< 


I 




1 





R{u)dT 



X' 



, + ^ \\R{u)\\^s + a,,-l. 



We will prove the following proposition in the next section. 
Proposition 1. For fixed s > and a < min(2s, |), we have 

\\R{u)\\xs + a,b-l < \\u\Ws,b. 

provided that < 6 — | zs sufficiently small. The same inequality holds for the restricted 
norms. 

Using Proposition [T], we have for < t < 5 < 1 (where [0,5] is the local existence 
interval) 

(2) \\u{t)-e^''^^^g\\Hs+a 



< 



\uiT)\\%.dT + \\u\\is,t < \\uf^s,t < \\g\\Hs 



^6 



We note that (see page 113-114 in [7j) the local existence time 6 depends on the norm 
of the data, 6 = S{\\g\\2) < 1, and by iteration we also have for any s > 0, 

\\u\\hs < Ce^l*l||^||H- =: T{t). 

Fix t large. For r < t, we have the bound 

Mr)\\Hs<Tir)<Tit). 

Using ([2]) repeatedly, we have 

\\uij5) - e^'^-uiu - mus.. < Mu - mwi. < ^(t)^ 
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for any j < t/6. Using this we obtain (with J = t/6) 

J 

\\u{J6) - e'^^^-^f\\Hs+a < J2 \\e'^^~'^^^''''u{j6) - e"^^-^^^^^^-u{{j - 

i=i 
J 

This finishes the proof of Theorem |2l 

4. Proof of Proposition [T] 
We start with the following elementary lemma (see, e.g., the Appendix of [13j). 

Lemma 1. // /3 > 7 > and /3 + 7 > 1, then 

E 7 rWl Tw ~ - - h), 

^-^ in — ki)P{n — kor 

where 

^ f 1, /3>1, 

[ {ky-^ /3 < 1. 

Using the definition of X^'^ norm, we have 

{kY'^''u{ki, Ti)u{k2, T2)u{k ~ ki + k2,T - Ti + T2) 2 



\\R{u)\\xs+a^b-i - / E (r-k^V-b 

Let 

f{k,T) = \u{k,T)\{kr{r-ky. 

It suffices to prove that 

/ E M{ki,k2,k,Ti,T2,T)f{ki,Ti)f{k2,T2)f{k-ki + k2,T-ri+T2) 



where 

M{ki,k2,k,Ti,T2,T) 



< 11/11^ = 11^115..^, 



{kY+^{ki)-'^{k2)-'{k-ki + k2y 



(r - A;2)i-^'(ri - klY{T2 - kj^ir - n + T2 - {k - ki + k2YY 
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By Cauchy-Schwarz in Ti,T2,ki, k2 variables, we estimate the norm above by 



■P(/ Yl M\h,k2,k,n,T2,T)')x 

I Yl /'(^l' T2)fik -k, + k2,T-n + T2] 



1 r 1 



Note that the norm above is equal to 11/^ * * /^ILi ri , which can be estimated by 
by Young's inequality. Therefore, it suffices to prove that the supremum above is finite. 
Using Lemma [1] in ri and T2 integrals, we obtain 

/ \- M^< \- {kr^+'^{k,)-'^k2)-'^{k-k,+k2)-'^ 

{ky'+^''{k,)~^'{k2y^'{k -ki + k2)~^' 



< 



sup 



The last line follows by the simple fact 

(3) (r — n){T — m) > {n — m). 

Since we have only the nonresonant terms, it suffices to estimate 

{k - A;i)2-2fe(A;i - fc2)2-2f' 
To estimate this sum we consider the cases \k — ki + kol > \k\, l/cil > \k\, and Ikol > \k\. 

I J. I Z( I I I J I -L I (-NJ I I J 1^1 rN^/ I I 

In the estimate we always take < 6 — | sufficiently small for any given s > and 
< a < min(2s, |). 

In the first case, using Lemma [1] we bound the sum by 



{k - kiY-^\ki - ^2)2-2?' ~ ^ (A; - A;i)2-26(fc^)2s+min(2-2fe,2s) ' 



In the case s > ^, we bound the sum by 

(A;)2'^log(l- 
(A;- A:i)2-2^>(A;i)2«+2 



^ (A;)2"log(l+(fci)) ^ . 2a-2+26+ < . 
/t- _ i--,\2-26/l. \2s+2-26 ~ \ ' ~ 



provided a < i. 

In the case < s < using Lemma [1] we bound the sum by 

ikf'^ ^ J ^^)2a+4fe-4«-2 < S < i. 



(A: - A;,)2-26(A:,)4.+l-26 ~ I ^^^2a+2.-2 I < s < 
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This is bounded in k provided that < a < min(2s, |). 

The second case is identical to the first case after renaming the variables: rii = k — ki + k2, 
n2 = k2. 

In the third case, after renaming the variables rii = ki, n2 = k — ki + k2, and using 
Lemma [1] we bound the sum by 

(^k - ni)2-26(A; - n2)2-2^ ~ ^ ' '^max{2s,2-2f,)l«^J\«/ 



ni,n2 



provided that s > and < a < min(2s, 1). This finishes the proof of Proposition [T] 



Remark 2. A similar smoothing property was established in [12] for the KdV equation. 
However, in that case, smoothing is not immediate within the context of the X^'^ norms as 
the following example shows [15]. Fix s G M. Consider u{n, r) = 5{n — M)x[-i,i]{t — n^), 
and v{n, r) = 6{n — l)x[_i^i](r — n^) with M 1. The inequality 

\\{uv)x\\yS+a,~l/2 < .,l/2||f||^s,l/2 
^KdV ^KdV ^KdV 

fails for any a > 0. 

The analogue of Theorem [2] was proved in [12] after transforming the KdV equation 
using differentiation by parts. However, the dimension statement in Theorem [1] is open for 
KdV since the corresponding result for the linear part (Airy equation) is not known. We 
will address this issue in future work. 

Remark 3. The gain of almost half a derivative in Theorem [2] can be inferred by taking 
the first Picard iteration for the cubic NLS and test the term 

Jo 

On the Fourier side, ignoring resonant terms, one sees that for g & L"^ the term 

r e-^"^^'+'''^'^'''+'''^''g{ki)g{k2)g{k3)dt' 

^ ^ g{ki)g{k2)g{k3) 2i(k^+k2)(k2+k,)t .^ 

is in but not necessarily better. 
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